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ABSTRACT 

We present a 63 ks FUSE observation of the Narrow-Line Seyfert 1 galaxy Arakelian 564. The spectrum is 
dominated by the strong emission in the O VI AA1032, 1038 resonance doublet. Strong, heavily saturated ab- 
sorption troughs due to Lyman series of Hydrogen, O VI and C III A977 at velocities near the systemic redshift 
of Arakelian 564 are also observed. We used the column densities of O VI and C III in conjunction with the 
published column densities of species observed in the UV and X-ray bands to derive constraints on the physical 
parameters of the absorber through photoionization modeling. The available data suggest that the UV and X-ray 
absorbers in Arakelian 564 are physically related, and possibly identical. The combination of constraints indicates 
that the absorber is characterized by a narrow range in total column density Nn and U, centered at logVe ~ 21 
and \ogU rs —1.5, and may be spatially extended along the line of sight. 

Subject headings: galaxies: active - galaxies: individual (Arakelian 564) - galaxies: nuclei - galaxies: Seyfert - 
FUV: Seyfert 



with z = 0.02467 and V = 14.6 mag (de Vaucouleurs et al. 1991), 
and L2-iokeV = 2.4 x 10 43 ergs s" 1 (Turner et al. 2001, here- 
after Paper I). It was the object of an intense multiwavelength 
monitoring campaign that included simultaneous observations 
from ASCA (2000 June 1 to July 6, Paper I; Pounds et al. 2001; 
Edelson et al. 2002), HST (2000 May 9 to July 8, Collier et 
al. 2001, Paper II; Crenshaw et al. 2002, Paper IV) and from 
many ground-based observatories as part of an AGN Watch 4 
project (1998 Nov to 2001 Jan, Shemmer et al. 2001, Pa- 
per III). Akn 564 has shown a strong associated UV absorber 
(Crenshaw et al. 1999, Paper II, Paper IV). There are indica- 
tions that it also possesses a warm X-ray absorber, as seen by 
the absorption lines of O VII and O VIII detected in a Chandra 
spectrum (Matsumoto, Leighly, & Marshall 2001) 5 . 

In this paper we present the results from a 63 ks FUSE obser- 
vation of Akn 564 obtained on 2001 June 29-30 UT, focusing 
in particular on the O VI intrinsic absorption; we investigate the 
physical properties of the UV and X-ray absorbing gas using the 
constraints on column densities obtained during the multiwave- 
length observations of this AGN. In §2 we present the data. In 
§3 we describe our analysis methods. In §4 we test the hypoth- 
esis that the Warm UV-X-ray absorber are one and the same 
through photoionization calculations. In §5 we discuss some 
implications of our investigation. Our results are summarized 
in §6. In a forthcoming paper (Romano et al., in preparation) 
we will analyze the intrinsic SED of Ark 564 and the properties 
of the gas responsible for the broad emission lines. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed Ark 564 with FUSE (Moos et al. 2000; Sah- 
now et al. 2000) for 63 ks starting on 2001 June 29 07:37:42 
UT. The observations, consisting of 24 separate exposures, 
were performed in photon address (time-tag) mode through the 
30" x 30" low-resolution (LWRS) aperture. During our 

1 Based on observations made with the NASA-CNES-CSA Far Ultraviolet Spectroscopic Explorer. FUSE is operated for NASA by the Johns Hopkins University 
under NASA contract NAS5-32985. 

2 Department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, OH 43210-1173; promano, mathur, pogge, peterson@astronomy.ohio- 
state.edu. 

3 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, MS 83 Cambridge, MA 02138; jkuraszkiewicz@cfa.harvard.edu. 

4 All publicly available data and complete references to published AGN Watch papers can be found at http://www.astronomy.ohio-state.edu/~agnwatch. 

5 http://www.pha.jhu.edu/groups/astro/workshop2001/papers/. 
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1. INTRODUCTION 

More than half of the Seyfert 1 population shows optical/UV 
intrinsic absorption associated with their active nucleus (Cren- 
shaw et al. 1999, and references therein). The strong UV ab- 
sorption lines, Lya, C IV, N V (and less frequently Si IV and 
Mg II) are found to be blueshifted, or at rest, with respect to the 
narrow emission lines, providing an important indication of the 
presence of a net radial outflow of the absorbing gas. A similar 
percentage also shows an associated ionized ("warm") X-ray 
absorber (George et al. 1998; Reynolds 1997) that is character- 
ized by high ionization, U = 0.1-10 (U = Q/Airr 2 n-nc, where Q 
is the number of ionizing photons) and high total Hydrogen col- 
umn density, N n = 10 21 -10 23 cm , which signature is typically 
the presence of O VII and O VIII edges. During the last decade 
evidence has been accumulated indicating that the same gas is 
responsible for the absorption in the UV and X-ray spectra of 
Seyfert Is (Mathur 1994; Mathur et al. 1994; Mathur, Elvis, & 
Wilkes 1995; Crenshaw et al. 1999; Kriss et al. 2000; Monier 
et al. 2001; Brotherton et al. 2002). Although it is not always 
possible to model the X-ray/UV absorber as a single zone (es- 
pecially when the complex UV absorption is resolved in multi- 
ple velocity components that are characterized by a large range 
of column densities and ionization) common characteristics of 
these absorbers have emerged, i.e., they are composed of high 
ionization, low density, high column density gas that is out- 
flowing and is located in or outside the broad emission-line re- 
gion (BELR). It is therefore worthwhile to investigate the nature 
of this nuclear component in active galactic nuclei (AGN) that 
represents an outflow (or wind) that can carry away a signifi- 
cant amount of kinetic energy at a mass-loss rate comparable to 
the accretion rate required to fuel the AGN (Mathur, Elvis, & 
Wilkes 1995). 

Arakelian 564 (Ark 564, IRAS 22403+2927, MGC +05-53- 
012) is a bright, nearby Narrow-Line Seyfert 1 (NLS1) galaxy, 
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FIG. 1. — FUSE spectrum of Akn 564, binned to a resolution of 
0.6A (100 pixels). In addition to the prominent emission lines from the 
O VI AA1032, 1038 resonance doublet, we suggest identifications for the main 
emission and absorption features. The strong airglow lines (H I Lyman series, 
O I A989, O I A1027, and He I A584 seen in second order at 11 67 A), and 
the C III A977 (which is scattered solar light in the SiC detector) have been 
removed. All other absorption features are due to Galactic or intergalactic ab- 
sorption (indicated with "IS"). The absorption feature at ~ 1080A is partially 
due to a gap between detectors. 

observation, a high-voltage anomaly occurred, and Detector 1 
did not collect useful data during the first 8 exposures. To 
best model the background for this observation, we first used 
the task ttag_combine provided with the FUSE calibration 
pipeline, CalFUSE (version 2.0.5) 6 , to combine the last 16 ex- 
posures for Detector 1 (total exposure of 42 ks), and all 24 ex- 
posures for Detector 2. We then processed the combined ex- 
posures with the standard pipeline and extracted spectra from 
both detectors. The flux scale for the final spectra is accurate 
to ±10%, while the wavelength scale is accurate to ±15 km 
s" 1 . As a result of the high-voltage anomaly and data screen- 
ing, the effective on-source times were 41 ks in Detector 1A, 
39 ks in Detector IB, and 58 ks in Detector 2A, and 62 ks in 
Detector 2B. Consequently, the SiCl A and SiC IB spectra were 
discarded from further analysis, the final signal-to-noise ratio 
(S/N) being < 1.5 even at 0.6A (100 pixels) resolution. We 
also discarded the LiFIB spectrum since it showed wavelength- 
dependent differences in flux of up to 30-50 % compared to the 
LiFIA probably due to the the "worm", which cannot be cor- 
rected for by the pipeline (Oegerle, Murphy, & Kriss 2000). 

The full FUSE spectrum was obtained by combining the 
spectra extracted from the SiC2A, LiF2B, LiFIA, SiC2B, 
and LiF2A segments, yielding a wavelength coverage of 91 6- 
1175 A. We then rebinned the full spectrum in a linear wave- 
length scale using 0.07 A bins (10 pixels, here on our full- 
resolution spectrum, with an effective resolution of 20 km s" 1 ), 

6 See http://fuse.pha.jhu.edu/analysis/calfuse.html. 



0.2 A bins (30 pixels, medium-resolution spectrum), and 0.6A 
bins (100 pixels, low-resolution spectrum). 

Figure 1 shows the low-resolution spectrum after we cosmet- 
ically removed the strong airglow lines (mainly H I Lyman se- 
ries, O I A989, O I A1027, and He I A584 seen in second order at 
1 167 A) and C III A977, which is scattered solar light in the SiC 
detector. The main spectral features are identified, the most 
prominent being the emission lines of the O VI AA1032, 1038 
resonance doublet. Strong absorption features due to Ly/3 and 
O VI AA1032, 1038 at velocities near the redshift of Ark 564 
are also observed. We detect absorption from Galactic ISM 
molecular Hydrogen, mainly H2 Lyman series absorption (see 
Shull et al. 2000, Sembach et al. 2000, and Savage et al. 2000), 
and atomic Galactic ISM lines, including O VI AA1032, 1038 
(Mathur et al. 2002, in preparation). No intrinsic Lyman edge 
is detected. 

3. DATA ANALYSIS 

Our goal was to determine the column densities of the ionic 
species we observed in our spectrum, combine this informa- 
tion with the column densities available in the literature for 
Ark 564, and derive constraints on the physical parameters of 
the absorber (total density and ionization) through photoion- 
ization modeling. HST/FOS spectra of Ark 564 (Crenshaw 
et al. 1999) show the presence of strong intrinsic absorption 
lines of Lya, N VAA1238.8, 1242.8, Si IVAA1393. 8, 1402.8, 
and C IVAA1548.2, 1550.8. Of these lines, which are resolved 
in STIS spectra into multiple components (Paper II, Paper 
IV), Lya, N V, and C IV are completely saturated. Figure 2 
shows the FUSE full-resolution (0.07 A, 10 pixels) spectrum of 
Akn 564, in the Ly/3/O VI wavelength region. Close examina- 
tion of the O VI troughs shows that the lines are heavily satu- 
rated, and their shape is mainly determined by partial covering 
effects (see §3.2 and Figure 3 below). Therefore, the absorption 
lines are not resolved into components at different velocities 
with respect to the systemic velocity, and we must treat each of 
the absorption troughs as a single absorption component, and 
we can only determine the velocity-averaged column densities 
of the observed species. We note, however, that O VIA1032 is 
not completely black. Indeed, analysis of the spectra obtained 
with different pulse height restrictions and from night-only data 
(we did not use the latter for this work, since the lower S/N did 
not allow a proper subtraction of Galactic molecular Hydro- 
gen) shows that scattered light is marginal in this observation, 
and that there is no filling in of the absorption troughs. It is 
also clear that the uncertainty in the measurement of O VI ab- 
sorption line parameters, hence O VI column density, is dom- 
inated by the uncertainty in the underlying emission-line pro- 
file. Additionally, there is contamination from absorption lines 
of Galactic molecular Hydrogen, with a column density log 
Af(H 2 ) > 16 (K. R. Sembach 2002, private communication); 
this is not surprising, given the substantial amount of neutral 
atomic hydrogen (Nh = 6.4 x 10 20 cm" 2 , Dickey & Lockman 
1990) along the line of sight toward Ark 564. Given these lim- 
itations, we proceeded as follows: we determined a power-law 
continuum underlying the Ly/3/O VI wavelength region using 
the low-resolution spectrum (S/N < 15 in the continuum); we 
then modeled the Ly/3 and O VI emission lines from the highl- 
and medium-resolution spectra (S/N < 10 in the emission lines 
for high-resolution), and used H2 templates to estimate the H2 
contribution to the O VI absorption troughs. This part of the 
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FIG. 2. — Full-resolution (0.07 A, ~ 10 pixels) spectrum of Akn 564, in the 
Ly/3/O VI wavelength region. The dashed vertical lines mark the rest-frame 
wavelengths of the O VI lines and the ticks mark Galactic absorption lines. 
Overalid are our five adopted models for the combined continuum and emis- 
sion lines, low-lying (Low), medium (Medl and Med2), and high-lying (High, 
that best follows the O VI A1038 peak), and a model with skewed Gaussian 
emission lines (Skew). All absorption lines are saturated, but while Ly/3 and 
O VI A1038 are black, O VI A1032 is not. The mean 1-cr errorbar on the spec- 
trum is also shown for reference. A color version of this plot is available in the 
electronic edition of this Journal. 

analysis was done using the I RAF 7 task specfit (Kriss 
1994) in the STSDAS package. Finally, we measured the ab- 
sorption line parameters of the normalized line profiles using 
the high-resolution spectrum, and determined the column den- 
sities with the apparent optical depth method, which we briefly 
describe below (§3.2). 



3.1. Intrinsic O VI Emission Models 

In addition to the power-law continuum (that we kept fixed 
relative to the fit of the low-resolution spectrum) our models 
for the adopted "continuum" under the absorption troughs in- 
cluded a pair of broad O VI emission lines (FWHM = 4000- 
5000 km s" 1 ), a pair of narrow O VI emission lines (FWHM 
= 1000-1 100 km s" 1 ), and a broad and narrow Ly/3 emission 
lines. All lines are taken to have Gaussian profiles. The in- 
tensity of the O VI doublet lines was fixed to the optically thin 
value 2:1 for both broad and narrow lines, while their wave- 
lengths were linked to the ratio of their laboratory values; the 
FWHM and wavelength of the broad Ly/3 line were linked to 
the ones of the O VI lines, though a small shift in wavelength 
was permitted, consistent with increasing blueshift with respect 
to the systemic velocity as the ionization increases. 

To assess the possible range of absorption-line parameters, 
we considered five different models for the emission lines, 
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FIG. 3. — Top: the full-resolution spectrum of Akn 564, in the Ly/3/O VI 
wavelength region, with a model for the continuum, the emission lines (Medl) 
and molecular Hydrogen absorption lines overlaid as a solid thicker curve 
(logA'(H2) ~ 10 16 ). The other panels show the normalized line profiles of 
the absorption system at z a ~ Ze- The dashed vertical lines mark the rest-frame 
wavelengths of the O VI lines. 

shown in Figure 2, namely, low-lying (Low), medium (Medl 
and Med2), and high-lying (High, that best follows the O VI 
A1038 peak, though it clearly over-predicts the O VI A1032 
peak). In all four cases the Gaussian profiles for the emis- 
sion lines were symmetrical. As observed in many AGNs (e.g. 
Marziani et al. 1996), the emission lines in all our models are 
blueshifted with respect to the systemic redshift z e = 0.02467 as 
derived from H I measurements (de Vaucouleurs et al. 1991) by 
390-1240 km s" 1 . Finally, we considered a model in which the 
broad lines have the least blueshift (100 km s" 1 ) with respect to 
the systemic redshift, and the narrow lines are at the systemic 
redshift; in this case, in order to model the profile, the narrow 
emission lines must be highly asymmetrical (skewness = 0.2). 
The motivation for considering the last model is the increas- 
ing evidence that the high-ionization emission lines in NLSls 
are broader and present an excess of flux in the blue with re- 
spect to the low-ionization lines (Laor et al. 1997b; Peterson et 
al. 2000; Mathur 2000; Leighly 2001, and references therein). 
With our choice of Low, Medl, Med2, and High emission line 
profiles the absorption is redshifted with respect to the emis- 
sion lines. This is highly unusual, though not unprecedented 
(Mathur, Elvis, & Wilkes 1999; Goodrich 2000). When model- 
ing the emission with the Skew profile, the absorption is in part 
blueshifted, in part redshifted with respect to the emission lines. 
Table 1 lists the O VI model emission-line parameters. Since no 
inflection points are clearly seen in the observed emission line 
profiles, our models may not have a direct physical interpreta- 
tion. Also, while we believe the true shape of the emission line 



7 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under 
cooperative agreement with the National Science Foundation. 
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FIG. 4. — Top: normalized line profiles (Medl) of the absorption system as 
a function of radial velocity relative to a systemic redshift Zs = 0.02467 (H I). 
The arrows show the position of the zero velocity with respect to the broad 
emission lines (BEL) and narrow emission lines (NEL). Middle: covering fac- 
tor as a function of radial velocity. Bottom: optical depth as a function of 
radial velocity. The solid lines refers to the O VI A 1032 profile, the dotted line 
to the the O VI A1038 profile. 

profile may be most realistically represented by the Medl and 
Med2 profiles, we consider the most extreme profiles (Low, 
High and Skew) to bracket the constraints on the values of col- 
umn density. 

Our spectrum also shows narrow H2 and H I absorption-lines 
from the inter-stellar medium (ISM) in the Ly/3/O VI wave- 
length region; in particular, the top panel of Figure 3 shows six 
of such absorption lines that lie outside the absorption troughs, 
which we fit as Gaussians with specf it (FWHM = 30km 
s" 1 ). We then used H2 and H I absorption-line templates (S. Mc- 
Candliss 2001, private communication; K. R. Sembach 2002, 
private communication) to predict the position of the lines in the 
absorption troughs and derive their intensity by scaling them 
to the two lines at ~ 1065 A. Incidentally, we note a 0.1 15 A 
shift between the predicted and observed wavelengths of H2 
absorption lines. The amount of this shift is constant along our 
spectrum, hence we interpreted it as a residual zero-point offset 
in the wavelength calibration of our spectrum; as reported in 
many observations 8 , these offsets can be as high as 0.25 A in 
the LWRS aperture. We matched templates and spectra accord- 
ingly. Figure 3 shows the line profiles normalized with respect 
to the combined continuum, emission-line, and H 2 profiles, for 
the five different O VI emission-line models. 

3.2. Intrinsic O VI Absorption Measurements 

As mentioned above, the O VI lines are so heavily saturated 
that we must treat each of the absorption troughs as a single ab- 
sorption component (as opposed to the many components ob- 

8 http://fuse.pha.jhu.edu/analysis/calfuse_wpl.html. 



served in Si IV and Si IIIA1206.5; Paper IV). This assumption is 
only strictly valid if the physical conditions are approximately 
constant along the line profile, i.e., as a function of radial veloc- 
ity. We consired the possibility of partial covering of the lines 
and used the apparent optical depth method to determine the 
column densities. Following Hamann et al. (1997), we calcu- 
lated the lower limit to the line-of-sight covering factor Q from 
the residual intensities I r in the troughs as a function of radial 
velocity, and the corresponding apparent optical depth r, 



Q>l~y, 
If) 



t > In 



(1) 



(2) 



where Iq is the assumed continuum intensity across the absorp- 
tion line (in the case of O VI, Iq is our combined power-law con- 
tinuum and the emission line models corrected for the Galac- 
tic Hydrogen absorption as described in §3.1). Column densi- 
ties are then obtained by integrating the apparent optical depth 
across the line profile using (e.g., Savage & Sembach 1991) 



JVjon > 



ne 2 f\ 



t(v) dv 



(3) 



where A and / are the laboratory wavelength and oscillator 
strength of the transition, respectively. 

Figure 4 shows the normalized line profiles, the covering fac- 
tor, and the optical depth as a function of radial velocity relative 
to the systemic redshift z e = 0.02467 for the Medl emission line 
profile. Predictably, Q w 1 for most of the O VI A1038 profile 
and is consistent with unity for O VI A 1032. The arrows show 
the position of the zero velocity with respect to the broad emis- 
sion lines (BEL) and narrow emission lines (NEL). As noted in 
§3.1, while most of the absorption is blueshifted with respect 
to the systemic redshift, the absorption troughs are completely 
or at least partially (Skew model) redshifted with respect to the 
broad and narrow emission lines. The absorption troughs also 
show the presence of gas which is redshifted with respect to 
the systemic velocity. This may indicate that the absorbing gas 
is undergoing net radial infall, as is the case for NGC 5548 
(Mathur, Elvis, & Wilkes 1999) and RX JO 134-42 (Goodrich 
2000). 

We used log /A = 2.137 for O VI A1032 and log/A = 1.836 
for O VI A1038 (Morton 1991). Table 2 reports the val- 
ues of O VI column densities for our five assumed emission 
line models; N yi = [2.31,2.65] x 10 15 cm" 2 and [5.28,5.96] x 
10 15 cm" 2 when measured from O VI A1032 and O VI A1038, 
respectively. The errors quoted in Table 2 are relative to the 
measurement of the integral of t in velocity space only. We es- 
timate that molecular Hydrogen lines contribute <~ 10 % to the 
flux in the absorption troughs. We adopt (5.7 ± 0.07) x 10 15 
cm" 2 , obtained averaging the values from Medl and Med2 
emission line models for the O VI A1038 line, as a conserva- 
tive lower limit on the O VI column density. 

3.3. C III column density 

We also determined the column density of C III from the C III 
A977 absorption trough with the apparent optical depth method 
described in §3.2 (log /A = 2.872, Morton 1991). The Ly 7 /C III 
wavelength region does not necessarily require the same degree 
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FIG. 5. — Comparison of the adopted SEDs for Ark 564 normalized to the 
absorption-corrected rest-frame flux at 2 keV. The long-dash line is the "table 
AGN" model in Cloudy; the dotted line (LZ) is the ionizing continuum de- 
scribed in Laor et al. (1997a) and Zheng et al. (1997); the solid line (SED1) is 
the Ark 564 SED one described in Romano et al. 2002 (in preparation) which 
is based on data corrected for reddening with a standard Galactic extinction 
curve and E(B — V) = 0.06; the short-dash line (SED2) is the Ark 564 ionizing 
continuum which is based on data corrected for reddening with both a stan- 
dard Galactic extinction curve and E(B-V) = 0.03 and the intrinsic extinction 
curve Crenshaw et al. (2002) derive for Ark 564 and E(B-V) = 0.14. The 
square points denote adopted points for SED1, the circles for SED2, while the 
horizontal solid lines show the full ranges where data were available. A color 
version of this plot is available in the electronic edition of this Journal. 

of complication in the emission line profile as the Ly/3/O VI 
wavelength region; however, though we performed fits and 
measurements separate from the Ly/J/O VI ones, for consis- 
tency we adopted the same, albeit with different flux normaliza- 
tion, emission model as Med2: one broad (FWHM = 4000km 
s" 1 ) and one narrow (FWHM = 1100 km s" 1 ) Gaussian emis- 
sion line for C III and one broad Ly7 Gaussian emission line; 
the FWHM and wavelength of the broad Ly7 line were linked 
to the ones of the broad C III line. We obtain a C III column 
density of (3.19 ± 0.05) x 10 14 cm" 2 (errors are relative to the 
measurement of the integral of r in velocity space only, while 
a 10% contribution is due to molecular Hydrogen absorption 
lines). 

3.4. Velocity Centroids 

Table 3 reports the values of the radial velocity centroids rel- 
ative to the systemic redshift of O VI, C III, Ly/3, and Ly7 ab- 
sorption lines, along with the measured column densities. We 
also show the results of Paper IV to emphasize the good agree- 
ment between the centroid velocity shifts (relative to systemic 
redshift) obtained for the H I Lyman series in the FUSE and 
HST spectra. In Paper IV, it is also noted that while saturation 
is probably responsible for the discrepancy in the values of the 

9 http://www.pa.uky.edu/~gary/cloudy/. 



centroids in the different ions, the most saturated lines, i.e., the 
H I Lyman series, give us an estimate of the total coverage of the 
absorber. Paper IV reports a range [-420, +180] km s -1 forLya, 
and we obtain [-431,+174]km s" 1 for Ly/3, [-412,+147]km 
s" 1 forLy 7 , [-395, +177] km s" 1 for C III, and [-412,+130]km 
s _1 for O VI A1038. The less saturated O VI A1032 yields 
[-374, +142] km s . Our FUSE spectrum was obtained a year 
after the last of the HST/STIS spectra were taken, and we con- 
firm the finding of Paper IV that there are no changes in radial 
velocity coverage of the absorber. 

4. PHOTOIONIZATION MODELING 

It is common practice to use photoionization codes to pre- 
dict the fractional abundance of an element in a given ioniza- 
tion state, /i on , given an input continuum, density «, total col- 
umn density Nh, and ionization parameter U of the gas. The 
fractional abundance of an ion of an element X is related to 
its column density N\ oa and the abundance of its element Nx 
by N, on = NuNxNion, which provides a prediction of N on that 
can be tested against observations. We used Cloudy 9 (v94.00, 
Ferland 1996) to calculate /; on for the ionic species for which 
we measured column densities from the FUSE spectrum, and 
for the species with published column densities, which are re- 
ported for easy reference in Table 3. We considered a range of 
values of Nh for a range of input continua (described in detail in 
§4. 1), a total Hydrogen density of 10 5 cm -3 , and assumed solar 
abundances relative to Hydrogen. In the case of "table agn", we 
specified a grid of Nh and U values. For all other models, we 
normalized the SEDs with respect to the measured X-ray lumi- 
nosity in the absorption-corrected rest-frame 2-10keV energy 
range (L2-iokeV = 2.4 x 10 43 ergs s" 1 ; Paper I), and specified 
the radius of the cloud, thus obtaining U. We note that the use 
of the observed SEDs assumes that the absorbing gas sees the 
same ionizing continuum as the observer does. 

4. 1 . Input continua 

Figure 5 illustrates our choices of input continua for 

Cloudy. 

1. The Cloudy "table agn" continuum, which is the 
Mathews & Ferland (1987) continuum modified with a 
sub-millimeter break at 10 /im, so that the spectral index 
is changed from -1 to -5 /2 (specific flux F v oc v~ a ) for 
frequencies below the millimiter break. While "table 
agn" is unlikely to be a representative spectral energy 
distribution (SED) of Seyferts, we use this continuum 
for comparison with the literature. 

2. A combination of the SED described in Laor et al. 
(1997a) and Zheng et al. (1997) for radio-quiet objects 
(LZ in Figure 5); we have extended the original SED in 
Laor et al. (1997a) to cover the whole 10" 5 -7.354 x 10 6 
Ryd energy range as required by Cloudy. At the low 
energies, we defined a sub-millimiter break and at the 
high energies, a break at lOOkeV (with a spectral index 
of -5/3), analogous to the ones in "table agn". This 
"composite" continuum might be a typical AGN contin- 
uum. 

3. Observed SEDs: we used data obtained during the mul- 
tiwavelength monitoring campaign performed in 2000, 
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that included simultaneous observations of Ark 564 
from ASCA (Paper I), HST (Paper II), and from many 
ground-based observatories (Paper III). In addition, 
IRAS measurements (Moshir et al. 1990) and our FUSE 
observations have been used. While the full extent of 
the data is used to create a quasi-simultaneous SED (Ro- 
mano et al. in preparation), our adopted continuum for 
Cloudy only consists of selected points (also shown in 
Figure 5). All data have been corrected for redshift. 

Special care has been paid in correcting the data for red- 
dening, given the indications (Paper IV) that strong in- 
trinsic neutral absorption is present in Ark 564 in excess 
of the Galactic absorption. Therefore, we corrected the 
spectra for reddening in 2 different ways: 

(a) Using a standard Galactic extinction curve with 
E(B-V) = 0.06 mag (Schlegel, Finkbeiner, & 
Davis 1998) (SED1 in Figure 5); 

(b) Using a standard Galactic extinction curve with 
E(B-V) = 0.03 mag plus the intrinsic extinction 
curve Crenshaw et al. (2002) (Paper IV) derive 
for Ark 564 and E(B-V) = 0.14 mag (SED2 in 
Figure 5). In the FUSE band we extrapolated 
the extinction correction linearly from the one 
relative to the HST band. We also note that to 
match our FUSE spectrum and HST spectrum in 
the overlapping region, we had to scale the FUSE 
fluxes by 0.75. This is not inconsistent with a 
combination of effects such as flux intercalibra- 
tion uncertainties and, most importantly, source 
flux variability. 

In the X-ray, we used continuum points from the power- 
law fit (photon index V = a+l = 2.538) and added a 
black body component of temperature r=1.8xl0 6 K 
and luminosity Lbb = 2.48 x 10 38 ergs s" 1 , as derived 
from fits to the mean ASCA spectrum. The square points 
in Figure 5 denote the adopted points for SED1, the cir- 
cles for SED2, while the horizontal solid lines show the 
full ranges where data were available. The full details 
of the observed SED as well as the analysis of the prop- 
erties of the BELR gas will be presented in Romano et 
al. (in preparation). 

4.2. Physical conditions of the UV/X-ray absorber 

Following Arav et al. (2001) we constrain the characteristics 
of the absorber by plotting curves of constant A^n on the log U - 
logA^H plane. In this plane, for each constant Ni on curve, lower 
limits on the column densities, derived from apparent optical 
depth line fitting methods, exclude the area below it, while up- 
per limits, derived from non-detections, exclude the area above 
it. Figure 6 shows the A^ on constraints (see Table 3) for the 
"table AGN" input continuum, and solar abundances. Lower 
limits are shown as solid lines, upper limits as dashed lines, de- 
tections as dotted lines. The combination of constraints given 
by the column densities suggests that the absorber in Ark 564 
is characterized by a narrow range in A^ on and U, i.e., logt/ = 
[-1.74,-0.74] and logJV H = [19.90,21.89]. We note the con- 
sistency of all constraints without departure from solar abun- 
dances. Analogously, Figure 7 shows the Ni on constraints for 
our SED 1 input continuum and indicates log C/ = [—1.99, — 1.31] 




log ([/) 

FIG. 6. — Photoionization curves at constant ionic column density on the 
plane of total hydrogen column density Nu vs. ionization parameter U. The 
shape of the ionizing radiation is defined by the "table AGN" model in Cloudy, 
and the abundances are solar. Lower limits are shown as solid lines, upper lim- 
its as dashed lines, detections as dotted lines. For clarity, curves relative to 
ions of the same element have been drawn in the same color: C in red, O in 
light blue, Si in dark blue. The shaded region corresponds to the locus on the 
log CMog N}i space where all conditions are met (see §4.2), that is, where 
log U = [-1.74.-0.74], logAh = [19.90,21.89]. A color version of this plot is 
available in the electronic edition of this Journal. 




log (U) 

FIG. 7. — Same as Fig. 6, but with a ionizing continuum with no redden- 
ing intrinsic to Ark 564 (SED1 in Figure 5). Solar abundances are assumed. 
Most constraints are met in \ogU = [-1.99,-1.31], logjV H = [19.99,21.09] (see 
§4.2). 
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FIG. 8. — Same as Fig. 6, but with a ionizing continuum with Galactic 
and intrinsic reddening (SED2 in Figure 5). Solar abundances are assumed. 
Most constraints are met in log U = [-1.86,-1.02], logW H = [19.95,21.27] (see 
§4.2). The top axis indicates the distance of the absorbing gas from the con- 
tinuum source. 

and log/V H = [19.99,21.09], Figure 8 (SED2 continuum) indi- 
cates log!/ = [-1.86,-1.02] and log/V H = [19.95,21.27], while 
Figure 9 (LZ continuum) indicates \ogU = [-1.97,-1.54] and 
logA^H = [20.00,20.79]. Thus, depending on the input contin- 
uum, there are small but significant differences in the derived 
properties of the absorber. 

5. DISCUSSION 

The UV absorber in Ark 564 is in a general state of out- 
flow with respect to the systemic redshift (see Table 3). A very 
good agreement is found between the values of the velocity cen- 
troids we derive for the species observed in the FUSE spectrum 
and the ones derived for the HST/STIS spectrum (Paper IV); 
therefore, we adopt as the best estimate of the net radial ve- 
locity of the UV absorber the value obtained in Paper IV for 
Si III and Si IV, the least saturated lines: V out = -194 ±5 km 
s" 1 . The absorption troughs also show the presence of gas 
which is redshifted with respect to the systemic velocity. This 
can be explained in part as a saturation effect, as is shown in 
Figure 3 of Paper IV. Alternatively, a model with more than 
one kinematic component is required to explain the observed 
absorption troughs (i.e. Elvis 2000); in this scenario, in addi- 
tion to the blueshifted absorption from an outflowing wind, we 
would be observing redshifted absorption from infalling ma- 
terial, such as an accretion flow. In addition to the contin- 
uum source, the absorbing gas must cover a substantial portion 
of the BELR, since the absorption troughs are much deeper 
than the continuum level. Assuming the identity of the UV 
and X-ray absorbing gas, we have used the column densities 
of the observed species to constrain the physical conditions of 
the absorber. For the most realistic SED (SED2), we obtained 




log ([/) 



FIG. 9. — Same as Fig. 6, but with a ionizing continuum described in 
Laor et al. (1997a) and Zheng et al. (1997) (ZL in Figure 5). Solar abun- 
dances are assumed. Most constraints are met in log U = [-1.97,-1.54] and 
\ogN H = [20.00, 20.79] (see §4.2). 

log/V H = [19.95,21.27] and \ogU = [-1.86,-1.02]. These con- 
straints can be used to determine the size of the absorber, its dis- 
tance from the central continuum source, and the mass outflow 
rate. For the following order-of-magnitude arguments we adopt 
the mean values U = 0.0363 and N u = 4.07 x 10 20 cm -2 . The 
size of the absorber is derived from the total column density, 
Tabs = 4.07 x 10 20 «g cm. For SED2 the number of ionizing 
photons is Q = 6.68 x 10 55 s" 1 , so the distance from the con- 
tinuum source is R^s = 7.00 x I0 22 n^ 2 cm. Using the lower 
limit on R R \, S > 95 pc found in Paper IV, this would indicate 
a total density zih > 5.70 x 10 4 cm" 3 . Assuming uniform den- 
sity, and considering that r a b s <C ^ a bs, the mass of the outflow- 
ing gas is M abs =2.11 x 10 10 /"h ^0, where /is the covering 
factor, i.e., the fraction of the sky covered by the absorber as 
seen at the central source. The mass outflow rate is then M a b s = 
Mabs V ol ,t A a bs = 3.17 x 10 4 / Mq yr" 1 and the outflow carries 
out a rate of kinetic energy M a b s V 2 ut /2 = 3.76 x 10 44 /ergs 
s . To power Ark 564 at the observed luminosity (Lboi = 
10 x L2-iokeV = 2.4 x 10 44 ergs s -1 ) at an efficiency r\ = 0.1, an 
accretion rate M acc = 1.8 x 10" 3 (L44/77) =4.3 x 10" 2 M Q yr _1 is 
required (L44 is the bolometric luminosity in units of 10 44 ergs 
s" 1 ). 

The outflow carries out a kinetic luminosity about one or- 
der of magnitude smaller than the observed radiative lumi- 
nosity of the source. However, the mass outflow rate is un- 
comfortably large unless the covering factor is very small. If 
M d bs ^5 M xc , then it implies / < 10~ 6 . Alternatively, our as- 
sumption r a bs <C Rabs might not be valid. The absorber might 
be an extended, low density region. The assumption of a uni- 
form density gas may not be strictly valid and the ionization 
parameter and density that we deduced should only be consid- 
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ered as "average" values. Recent Chandra observations have 
found extended warm gas in some AGN (Sako et al. 2000) with 
physical characteristics similar to that of a warm absorber, but 
seen in emission. So it is quite likely that the warm absorber in 
Ark 564 is also spatially extended along the line of sight. 

In Paper IV the UV absorber was modeled as a single zone 
with quasi-solar abundances (Carbon depletion being the main 
departure) and the best-fit values of log C/ = —1 .48 and logVe = 
21 .21 are consistent with our limits. Crenshaw et al. (2002, Pa- 
per IV) over-predicted Carbon and Oxygen column densities: 
Afcm = 5.2 x 10 15 cm" 2 , V 0V i = 2.4 x 10 17 cm" 2 (cf. our mea- 
surements: Wcm = 3.2 x 10 14 cm -2 and AW = 5.7 x 10 15 cm -2 ). 
These predictions, however, are consistent with the upper-end 
values of our range of parameter space. We note that our mod- 
eling did not require Carbon to be depleted. 

Finally, we can compare our solutions of logA^ = 21 and 
\ogU = -1.5, with the preliminary results of Matsumoto, 
Leighly, & Marshall (2001) based on analysis of a 50 ks Chan- 
dra observation of Ark 564. Their curve of growth analy- 
sis on the absorption lines indicates that A/ovn = 3.2 x 10 17 
cm -2 , Novm = 1 x 10 18 cm -2 , A^eix = 3.2 x 10 17 cm -2 , and 
A , NeX= 1 x 10 17 cm" 2 , suggestive oflogA f H = 21 andlog£= 1.6- 
2. While there is agreement between the values of logA^ and 
No vn, the column densities they measure for O VII do not agree 
with the ones derived in Paper IV from the upper limits on 
the bound-free optical depths in the ASCA spectrum (Paper I). 
Given the high Afovni we would expect an edge would be ob- 
servable in the Chandra spectrum. 

6. SUMMARY 

We have presented a 63 ks FUSE observation of the NLS 1 
galaxy Ark 564. The observed spectrum is dominated by the 
O VI AA1032, 1038 emission lines. As observed in many AGNs 
(e.g. Marziani et al. 1996), the emission lines in all our mod- 
els are blueshifted (or at rest, as in the case of our Skew blue- 
asymmetric model) with respect to the systemic redshift by 
100-1240 km s" 1 . Blue-asymmetric UV emission line profiles 
may be a characteristic of NLS1 galaxies. 

We concentrated on the analysis of the strong and heavily 



saturated absorption troughs due to Lyman series, O VI and 
C III A977, which are observed at velocities near the systemic 
redshift of Ark 564. In a forthcoming paper (Romano et al., 
in preparation) we will analyze the intrinsic SED of Ark 564 
and the properties of the BELR gas. Using the apparent optical 
depth method, we have determined that the column density of 
O VI is a few 10 15 cm" 2 , and N C m = 3.2 x 10 14 cm" 2 . We used 
these values in conjunction with the published column densities 
of species observed in the UV and X-ray spectra of this object 
to derive constraints on the physical parameters of the UV/X- 
ray absorbing gas through photoionization modeling. The com- 
bination of constraints, assuming the most realistic SED, indi- 
cates that the absorber is characterized by a narrow range of 
density and ionization parameter, logA^ = [19.95,21.27] and 
logf/ = [-1.86,-1.02]. 

There is excellent agreement between the kinematic proper- 
ties of the UV absorber emerging from the combined analysis 
of the FUSE and HS77STIS spectra, i.e. distribution of gas in 
radial velocity (as derived from the extent of the absorption 
troughs) and net radial velocity (as derived from the velocity 
centroids). The UV/X-ray absorber in Ark 564 is in outflow 
with respect to the systemic redshift with a radial velocity of 
a few hundred km s" 1 , and it is likely spatially extended along 
the line of sight. The absorption troughs also show the presence 
of gas which is redshifted with respect to the systemic veloc- 
ity. This may indicate that a component in the absorbing gas 
is undergoing net radial infall. This is highly unusual, though 
not unprecedented (Mathur, Elvis, & Wilkes 1999; Goodrich 
2000). 
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Table 1 



O VI Model Emission-Line Parameters 



Model 


Line 


A a 


Flux 3 


FWHM 


AV b 






(A) 


(lO-^ergs" 1 cm" 2 ) 


(km s" 1 ) 


(km s" 1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Low 


BEL O VI A 1032 


1053. 


1.00 


5000. 


-1240 




BEL O VI A1038 


1058.82 


0.50 


5000. 


-1240 




NELO VI A 1032 


1056. 


0.90 


1100. 


-390 




NELO VI A1038 


1061.83 


0.45 


1100. 


-390 


Medl 


BEL O VI A 1032 


1055. 


1.00 


4000. 


-670 




BEL O VI A1038 


1060.83 


0.50 


4000. 


-670 




NELO VI A 1032 


1056. 


1.00 


1100. 


-390 




NELO VI A1038 


1061.83 


0.50 


1100. 


-390 


Med2 


BEL O VI A 1032 


1056. 


1.00 


4000. 


-390 




BEL O VI A1038 


1061.83 


0.50 


4000. 


-390 




NELO VI A 1032 


1056.3 


1.00 


1100. 


-300 




NELO VI A1038 


1062.13 


0.50 


1100. 


-300 


High 


BEL O VI A 1032 


1055. 


1.00 


4000. 


-670 




BEL O VI A1038 


1060.83 


0.50 


4000. 


-670 




NELO VI A 1032 


1056. 


1.25 


1000. 


-390 




NELO VI A1038 


1061.83 


0.63 


1000. 


-390 


Skew 


BEL O VI A 1032 


1057. 


1.25 


2000. 


-100 




BEL O VI A1038 


1062.84 


0.63 


2000. 


-100 




NELO VI A 1032 


1057.37 


0.60 


1000. 


0. 




NELO VI A1038 


1063.21 


0.30 


1000. 


0. 



a Observed values. 

b Velocities are relative to the systemic redshift z e = 0.02467 (H I measurements, de 
Vaucouleurs et al. 1991). The shift toward longer wavelengths of 0.115 A to match the 
Galactic molecular Hydrogen templates is not included. 
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Table 2 



O VI Column Densities from Intrinsic Absorption 



Line 


Low Model 


Medl Model 


Med2 Model 


High Model 


Skew Model 




(10 15 cm" 2 ) 


(10 15 cm" 2 ) 


(10 15 cm" 2 ) 


(10 15 cm" 2 ) 


(10 15 cm" 2 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


OVI A1032 


2.31 ±0.03 


2.47 ± 0.05 


2.54 ± 0.03 


2.65 ± 0.03 


2.80 ± 0.03 


OVI A1038 


5.28 ± 0.05 


5.60 ± 0.05 


5.77 ± 0.05 


5.96 ± 0.05 


5.64 ± 0.05 



Note. — The errors quoted are relative to the measurement of the integral of t in velocity 
space only. We estimate that molecular Hydrogen lines contributes ~ 10% to the flux in the 
absorption troughs. 



Table 3 

Column Densities from Intrinsic Absorption in Arakelian 564 



Ion 


Wavelength/ 


Lower Limit 


Detection 


Upper Limit 


Velocity 


References 




Energy 


(cm -2 ) 


(cm" 2 ) 


(cm" 2 ) 


(km s" 1 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


HI 


973 A 








-108 


1 


cm 


977.0 A 


3.2 x 10 14 






-153 


1 


HI 


1025 A 








-111 


1 


ovi a 


1031.9,1037.6A 


5.7 x 10 15 






-79/- 116 


1 


Si III 


1206.5 A 




2.6 x 10 13 




-190 


2 


HI 


1216A 


1.4 x 10 15 






-106 


2 


N V 


1238.8,1242.8A 


3.1 x 10 15 






-152 


2 


Sill 


1260.4 A 






7.4 x 10 13 




2 


CII 


1334.5A 






5.4 x 10 13 




2 


Si IV 


1393.8,1402.8A 




1.6 x 10 14 




-197 


2 


CIV 


1548.2,1550.8A 


2.5 x 10 15 






-130 


2 


ovn b 


0.74 keV 






2.2 x 10 17 




2 


o vm b 


0.87keV 






1.1 x 10 16 




2 



Note. — Col. (1): Ion. Col. (2): Wavelength of absorption lines or energy of absorption edges. Col. 

(3) : Lower limit on column densities as derived from line fitting/apparent optical depth methods. Col. 

(4) : Column density values derived using the multiplet method of Hamann et al. (1997). Col. (5): Upper 
limit on column densities as derived from non-detection of the line. Col. (6): Velocity centroids relative 
to the systemic redshift z e = 0.02467. Col. (7): References for column densities and velocity centroids. 

a Average of values obtained from Medl and Med2 emission line models. 
b Absorption edge. 
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